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A  Modeling  Study  of  the  Mesocyclone 


1.  INTRODUCTION 

A  well-known  characteiistic  of  supeicell-type  stonns  is  a  deep,  wide  region  ( >  3  km)  of  persistent  ( minutes) 
rotation  (shear  >  0.(X)S/s)  known  as  the  mesocyclone.  Strams  exhibiting  this  degree  of  organized  rotation  are  extremely 
dangerous.  Almost  all  produce  severe  weather  of  one  kind  or  another.  It  is  fortuitous  that  the  mesocyclone  produces  a 
distinct  signature  in  the  single-Dopplo'  velocity  field  that  is  not  only  easy  to  recognize  by  eye,  but  is  generally 
straightforward  to  detect  dirough  automated  means.  The  single-Dopplor  mesocyclone  signature  is  now  considered  to  be 
a  most  valuable  warning  criterion  for  severe  weather. 

Doppler  radar  only  senses  the  component  of  motion  that  is  directed  along  the  radar  beam.  This  limited  perspective 
means  that  a  rotating  mass,  larger  in  width  than  the  radar  beam,  will  be  sensed  directly  only  along  a  cross  section.  With 
only  one  dimension  of  the  rotating  feature  well  resolved,  the  total  shape  and  velocity  distribution  of  the  rotating  area  is 
not  easily  obtained  For  convenience,  and  for  lack  of  a  better  model,  the  mesocyclone  is  generally  assumed  to  be 
circular,  with  the  flow  varying  according  to  a  Ranldne-combined  vortex  (RCV)  as  first  proposed  by  Lhennitte  (1964), 
herein  referred  to  as  the  circular  mesocyclone  model 

The  assumption  of  circular  flow  is  a  reasonable  first-order  ^^oximation,  and  is  suppmted  by  numerous 
observations.  Circularity  is,  in  fact,  a  cornerstone  of  existing  detection  and  evaluation  techniques.  Donaldson  (1970) 
proposed  criteria^  for  positive  single-Dc^pler  identification  of  a  mesocyclone  based  on  a  circularity  assumption.  These 
criteria,  having  proven  effective  during  iterations,  were  built  in  one  form  or  another  into  automated  mesocyclone 
detection  algoriduns  (for  example,  Vfieler,^  and  Zmic’  et  al.,^ ).  Circularity  is  assumed  in  automated  techniques  for  hail 
prediction  (for  example,  Witt  and  Nelson,^)  and  tornado  prediction  (Desrochers  and  Donaldson,^).  The  success  of  severe 
weather  evaluates  that  use  mesocyclone  related  parameters  depends  largely  on  die  validity  of  the  model  used  to  extract 
the  mesocyclone  information. 

Mesocyclones  are  not  necessarily  circular.  There  is  evidence  to  support  this  in  multiple-Doppler  studies  (for 
example,  Brandes,^),  where  the  mesocyclone  sometimes  rqipears  elongated.  Mesocyclone  signatures  from 
single-Doppler  observations  often  differ  in  detail  from  the  ideal  proposed  by  the  circular  model.  Some  of  these 
irregularities  can  be  explained  by  embedded  vortices  within  a  general  mesocyclone  flow,  as  modeled  nicely  by  Brown 
and  Wood^.  Other  irregularities  are  more  difficult  to  explain,  except  to  say  they  must  be  associated  widi  a  distorted  or 
non-uniform  flow  field.  The  results  of  such  distortions  is  die  tc^ic  of  this  investigation. 
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As  a  starting  point,  it  is  useful  to  examine  the  circular-model  kinematics  to  see  how  well  they  agree  with  the 
three-dimensional  (3-D)  flow  field  of  a  superceU-type  updraft,  from  dure,  alternative  models  of  the  mesocyclone  will 
be  examined.  Specifically,  we  will  examine  an  elliptical  flow  field  to  show  how  it  can  better  explain  some  of  the 
irregularities  observed  in  the  single-Doppler  flow. 

2.  THE  CLASSICAL  MESOCYCLONE  MODEL 

The  velocity  distribution  of  the  circular  mesocyclone  model  is  shown  in  Figure  la.  This  model  consists  of  two  flow 
regimes:  an  inner  or  “core”  region,  and  an  outer  or  “potential  flow”.  Flow  within  the  core  varies  linearly  with  distance 
from  the  mesocyclone  center  and  is  described  by 

v/r  =  Ci,  (1) 

where  (v)  is  the  velocity  at  any  radius  (r)  inside  the  ewe,  and  Ci  =  V  /  R,  a  constant  related  by  die  peak  velocity  (V) 
found  at  the  mesocyclone  cote  radius  (R)  that  separates  the  solid  body  and  potendal  flows.  Conversely,  flow  within  the 
potential  flow  region  is  given  by 

vr  =  C2=VR.  (2) 

When  radar  resolution  permits,  (as  it  generally  does  for  flows  the  size  of  mesocyclones)  the  mesocyclone  is  seen  as 
a  velocity  couplet,  with  the  flow  defined  well  at  only  the  two  points  where  the  radar-relative  incoming  and  outgoing 
velocities  are  greatest.  An  example  of  a  couplet  for  pure  rotation  is  shown  in  Figure  lb,  where  the  radar  viewpoint  is 
from  the  bottom  of  the  page.  For  rotation,  the  model  describes  a  flow  that  rotates  as  a  solid  body  inside  the  core  region. 
Pure  rotation  may  exist  at  one  or  mote  levels  between  2  and  7  km  AGL. 

At  most  elevations  there  is  a  component  of  divergence  associated  with  the  mesocyclonic  flow.  This  results  in  a 
clockwise  rotation  of  the  couplet  in  Figure  lb  for  convergence  and  counterclockwise  for  divergence.  The  orientation 
angle  of  the  velocity  peaks  relative  to  the  radar  beam,  is  0®  for  pure  rotation,  -90®  for  pure  convergence,  and  90®  for 
pure  divergence.  At  the  mesocyclone  core  radius,  conqronents  of  velocity  associated  with  rotation  (Vr)  and  divergence 
(Vd)  are  given  by  Vo  cos  (<t>)  arid  Vo  sin(<^),  respectively,  where  Vo  is  the  average  of  the  radar  observed  velocity  couplet 
velocity  maxima  Figure  2  shows  a  mesocyclone  with  ^  =  -45°,  where  the  component  of  velocity  due  to  inflow  is  equal 
to  the  component  associated  with  rotation. 

Pure  convergence  may  exist  at  low  levels  of  a  storm,  but  genraally  not  as  part  of  a  circular,  mesocyclonic  flow.  At 
the  storm  top,  however,  signatures  of  pure  divergence  in  an  apparently  circular  flow  are  common. 

Vorticity  and  divergence  fields  are  simple  and  straightforward  to  evaluate  for  a  circular  mesocyclone  flow.  The 
circular  model  equations  describe  a  flow  field  in  which  magnitudes  of  velocity  are  a  function  only  of  range  from  the 
vortex  center  and  not  a  function  of  direction.  Consider  first  a  purely  rotational  flow  where  the  velocity  varies  according 
to  Eqs.(l)  and  (2)  and  the  total  component  of  velocity  is  normal  to  the  mesocyclone  radius  arm.  With  the  mesocyclone 
examined  in  a  Cartesian  corudinate  system  defined  by  directions  x  and  y,  and  velocity  components  Ux  and  uy,  then,  in 
the  mesocyclone  core,  duy/dx  is  positive  at  any  given  point  for  cyclonic  rotation,  and  equal  in  magnitude  but  opposite  in 
sign  to  9ux/9y.  Therefore,  within  the  core  region. 


Vorticity  ( VxV  )  =2  duy/dx  =  2  Vo  cos(<t>)/R  =  2  Vr/R. 

Vorticity  is  constant  within  the  core  region  and  equal  to  twice  the  peak  mesocyclone  rotational  velocity  divided  by  the 
core  radius. 
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Figure  1.  Velocity  Distribution  for  a  Circular,  Axisynunetric,  Rankine-Combined  Vortex  (RCV)  Mesocyclone  Flow. 

a)  Velocity  magnitude  as  a  function  of  distance  from  the  mesocyclone  center.  The  mesocyclone  radius,  called 
the  core  radius  (R),  corresponds  to  the  locations  of  peak  velocity  directed  towards  (-)  and  away  from  (+)  the 
radar  (±V).  b)  Them-edcal  mesocyclone  appearance  for  a  purely  rotational  flow  as  seen  by  single  Dof^ler  ra¬ 
dar,  for  a  radar  located  in  the  direction  towards  the  bottom  of  the  page.  Contour  levels  are  as  indicated  except 
that  radial  velocities  within  7.S  percent  of  maximum  or  minimum  are  indicated  by  the  solid  black  contour. 
The  circle  drawn  indicates  the  mesocyclone  ring  of  maximum  velocity  that  defims  the  c<»e.  Pure  rotation  is 
indicated  by  inbound  and  outbound  velocity  peaks  that  are  at  the  same  range  from  the  radar.  This  is  defined 
by  a  velocity  peaks’  orientation  angle  ((|)}  of  0°. 
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Figure  2.  Theoretical  Appearance,  as  Seen  by  Single-IX^pler  Radar,  of  a  Circular  Mosocyclone  Containing  Equal 
Magnitudes  of  Vcxticity  and  Convergence.  Contour  levels  are  as  indicated  in  Figure  1.  The  relative  magni¬ 
tudes  of  vorticity  to  convergence  or  divergence  is  indicated  by  the  velocity  peaks’  orientation  angle  (^)  meas¬ 
ured  relative  to  the  azimuthal  direction.  Here,  <|>  equals  -45°. 
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In  the  region  of  potential  flow  for  rotation,  Ux  and  u .  vary  with  distance  from  die  vortex  center  accmding  to  Eq.  (2). 
The  rotational  velocity  decreases  with  increasing  distaiK^e  from  the  mesocyclone  center  and  duy/dx  =  dux/dy.  The 
derivatives  are  equal  in  IxMh  sign  and  magnitude  at  any  point,  the  rate  of  change  of  each  (tecieases  with  increasing 
distance,  and  the  vector  {^oduct  for  vordcity  is  everywhere  equal  to  zero  in  the  potential  flow  area. 

It  seems  unlikely  that  vorticity  should  actually  change  abrupUy  from  zero  to  a  constant  across  the  ring  of  maximum 
mesocyclone  velocity.  Regardless  of  the  validity  of  assuming  solid  body  rotation  in  the  core,  it  is  reasonable  to  assume 
that  a  transition  zone  of  increasing  cyclonic  voidcity  exists  in  the  so-called  “potential  flow  region”.  If  the  rate  of  change 
of  rotational  velocity  in  the  potential  flow  region  is  decreased  from  that  prescribed  by  the  circular  model  potential  flow 
Eq.  (2),  such  a  transition  is  provided  (not  shown). 

The  divergence  field  associated  with  inflow  or  outflow  is  analogous  to  the  vorticity  field  of  the  circular  model  for 
rotation,  that  is,  dux/dx  -  3uy/dy.  Therefore,  divergence  within  the  core  is  given  by 


Divergence  ( V-  =2  9ux/3x  =  2  Vo  sin((t>)/R  =  2  Vd/R. 

The  vorticity  and  divergence  field  for  the  circular  mesocyclone  example  of  Figure  2  is  shown  in  Figure  3.  Recall 
that  in  Figure  2,  where  <t>  =  -45°,  the  magnitudes  of  vorticity  and  convergence  inside  the  core  are  equal.  The  circular 
model  is  characterized  by  uniform  vorticity  and  divergence  (convergence  in  this  case)  associated  the  with  core  region. 
In  the  region  described  by  potential  flow,  vorticity  and  divergence  are  both  zero. 

In  surtunary,  the  circular  mesocyclone  model  is  characterized  by  uniform  vcaticity  and  divergence  in  the  core 
region.  The  relative  magnitudes  of  these  quantities  are  determined  by  the  orientation  of  the  velocity  peaks  with  respect 
to  the  radar  (^).  From  an  observer’s  point  of  view,  uniform  divergence  or  vorticity  in  a  single-Doppler  velocity  field  is 
indicated  by  straight  isodops  (velocity  contours),  as  seen  in  the  core  region  of  Figures  lb  and  2.  This  is  sometimes,  but 
not  always,  a  good  approximation  of  what  is  seen  in  nature.  Also,  from  Figures  lb  and  2,  it  is  obvious  that  the  zero  line 
(indicat^  by  the  non-shaded  region)  is  always  perpendicular  to  a  line  connecting  the  velocity  peaks  in  the  core.  There  is 
generally  some  deviation  from  this  in  single-Doppler  observations.  In  the  “potential  flow”  region,  the  circular  model  is 
perhaps  less  accurate  than  for  the  core.  The  impact  of  this  for  evaluating  mesocyclone  intensity  is  minimal  because 
those  evaluations  concern  only  the  core,  where  the  concentration  of  energy  is  greatest. 

3.  NON-CIRCULAR  ROTATION  AND  NON-SYMMETRICAL  FLOWS 


It  has  been  recognized  that  mesocyclones  do  not  always  assume  a  circular,  RCV-like  flow  pattern.  Wood  (1S>91), 
for  example,  examined  the  envelope  of  azimuthal  gradients  in  the  single-Doppler  field  and  compared  it  to  the  envelope 
predicted  by  a  circular  flow^.  His  work  documents  several  single-Doppler  cases  where  the  mesocyclone  flow  appears  to 
deviate  considerably  firom  the  circular  mesocyclone  model.  The  observed  deviation  is  sometimes  due  to  embedded 
vortices  that  confuse  the  issue.  At  other  times  the  vortex  appears  singular  but  distorted  in  shape. 

In  the  multiple-Doppler  study  by  Brandes^  of  the  1977  Del  City,  Oklahoma,  mesocyclone,  where  the  total 
horizontal  wind  field  is  more  completely  determined,  the  mesocyclone,  as  defined  by  a  vorticity  magnitude,  is 
elongated.  An  example  from  this  work  will  be  discussed  in  Section  4. 

Deviations  from  the  circular  mesocyclone  model  are  sometimes  obvious  in  even  single-Doppler  data.  Distorted 
flows  are  seen  most  often  at  low  levels  but  can  extend  to  middle  storm  levels  as  well.  It  is  not  clear  at  this  time  what 
distortions  exist  at  the  upper  reaches  of  a  storm.  The  high  elevation  angles  required  for  sampling  these  regions  result  in 
degraded  horizontal  resolution  of  the  data  and  a  stronger  contribution  by  the  vertical  velocities  to  the  radial  component. 
Examples  of  highly  distorted  low  level  flows  from  the  Del  City  mesocyclone  are  shown  in  Figures  4  and  5.  The  times 
select^  are  for  volume  scans  beginning  at  1826  and  1832.  At  1840  the  storm’s  first  tornado  occurred,  classified  as  F2 
on  the  Fujita  scale*®. 
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Figure  3.  The  Associated  Vorticity  and  Divergence  (convergence)  Fields  for  a  Circular  RCV  Flow  Associated  with  the 
Example  in  Figure  2. 
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DEL  CITY  DEL  CITY 


B-scan  Format  of  the  Dei  City,  OK,  Storm  of  20  May  1977  at  1827  CST-  Data  as  Collected  From  the  National  Severe  Storms  Laboratory  10  cm  Dop 
pier  Radar  Located  in  Norman,  OK.  Range  and  azimuth  scaling  are  adjusted  so  they  represent  approximately  equal  distances.  Contour  levels  are  as 
indicated  in  Figure  1.  a)  The  0.7°  elevation  scan,  b) The  1.8° elevation  scan. 


The  scaling  piesented  in  Figures  4  and  S  is  nearly  equal  in  the  azimuthal  and  radial  dimensions.  Because  of  the 
range-azimuth  di^lays,  there  is  some  distotti<»  to  the  data.  However,  the  mesocyclone  is  sufficiently  far  from  the  radar 
that  the  distortion  attributable  to  the  display  is  small  compared  to  the  natural  distmtion  in  the  mesocyclone  flow.  Thoe 
are  several  distinctive  qualities  to  the  examples.  First,  variations  in  the  contouring  are  somewhat  suggestive  of  the 
circular  mesocyclone  model;  that  is,  there  appear  to  be  inner  and  outer  flow  regimes.  Second,  in  contrast  to  the  circular 
model,  the  line  connecting  the  velocity  peaks  in  the  Del  City  example  is  not  always  perpendicular  to  the  zero  line  (see 
Figures  4b  and  5b  for  exanq>le).  Finally,  the  envelope  in  the  isodops  separating  the  inner  and  outer  flow  regions  a{^>ears 
to  be  approximately  elliptical  in  stupe.  Hgure  Sb  illustrates  this  elliptical  sh^. 

A  principal  problem  to  solve  in  modeling  the  ell^cal  flow  is  the  determination  of  reasonable  kinematics 
constraints.  As  a  starting  point,  the  mesocyclone  is  examined  frtmi  a  two-dimensional  (2-D)  perspective,  although  it  is  in 
reality,  a  doee-dimensional  feature.  Thoefore,  constraints  can  only  be  applied  to  the  rotational  aspect  of  the  flow.  For 
simplicity,  the  model  is  developed  whoe  vorticity  is  unifmm  inside  the  core.  The  affect  of  applying  fields  of  unifmm 
divergence  is  then  examined.  The  process  used  for  combining  the  features  is  simple;  separate  fields  for  rotation  and 
inflow  are  generated;  and  dien  they  are  vectorially  combined  to  generate  the  total  two-dimensirmal  flow  field.  This 
simple  ^roach  allows  easy  examination  of  important  aspects  of  the  model. 


3.1  Rotafion  in  an  Elliptical  Flow 


The  orientation  of  the  Del  City  “ellipse”  with  respect  to  the  radar  is  relatively  unchanged  with  time  (Ingures  4  and 
S).  This  means  that  we  are  not  observing  an  elliptical  pattern  rotating  about  an  axis,  but  rather  an  elliptical  flow.  Also, 
the  changing  position  in  space  of  the  velocity  peaks  suggests  an  evolving  vorticity  and  divergence  field.  It  will  be  shown 
here  that  the  general  char^teristics  of  the  sin^e-Doppler  velocity  field  can  be  simulated  with  a  simple  two-dimensional 
(2-D)  mesocyclone  model. 

Mesocyclone  rotation  is  modeled  by  an  ellipse  with  votices  (±a),  covertices  (±b),  (as  shown  in  Figure  6a),  and 
thickness.  Cross-sectional  areas  along  a  and  b  through  which  the  flow  travels  are  given  by  A«  s  a  Ah  and  Ab  =  b  Ah  . 
Constant  mass  flux  is  assumed  for  the  rotational  flow  such  that 

pAaVa  =  pAbVb.  (5) 

where  Va  and  Vb  are  the  average  velocities  through  die  cross  sectional  areas,  and  p  is  air  density,  held  constant 
throughout  the  horizontal  flow.  For  a  rotational  flow  to  satisfy  Eq.  (5),  air  traveling  across  the  ellipse  major  axis  (a)  must 
experience  an  increase  in  angular  speed  (co)  to  pass  through  the  minor  axis  (b),  so  the  angular  speed  at  a  is  not  the  same 
as  b;  ooa  ^  (Ob .  If  0)  is  held  constant  along ^y  given  radius  arm  of  die  ellipse,  akin  to  the  inner  flow  of  the  circular 
mesocyclone  model,  then  (V  =  coa  a/2)  and  (Vb  =  (Ob  b/2).  Thoefore, 

pa2Ah<0a/2  =  pb2Ah<0b/2.  (6) 

The  peak  rotational  velocities  along  a  and  b  are  given  by  Va  and  Vb.  Therefore, 
coa  =  Va/a,  (Ob  =  Vb/b ,  and 


aVa  =  bVb.  (7) 

fri  a  purely  rotational  flow,  the  wind  is  tangent  to  the  ellipse  surface,  and  the  magnitude  of  the  rotational  velocity  is 
just  the  conqxment  of  the  velocity  that  is  normal  to  the  ellipse  radial  arm.  Along  axes  a  and  b  the  flow  tangent  to  the 
ellipse  is  also  normal  to  the  radius  arm,  but  for  all  radii  in  between,  the  rotational  velocity  is  given  by  the  velocity  along 
the  ellipse  surface  multiplied  by  sin(a),  where  a  is  die  obtuse  angle  between  the  radius  arm  of  the  ellipse  and  the 
tangent  to  the  ellipse  (Figure  d>).  Hie  ang^e  a  is  given  by  (90°  ^  -  P).  ^  is  the  ellipse  relative  angle  that  varies 
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between  0°  and  90°  with  0°  towards  the  ellipse  “top”  along  the  mimv  axis.  ^  is  the  arctangent  of  the  slope  along  the 
ellipse. 


Mth  an  ellipse  defined  by  (x^/a^  +  y^np-  =  1),  the  slope  along  the  curve  is  given  by 

dy _ -bx  (8) 

dx  a2(l 


If  die  coordinate  system  is  defined  with  \  measured  positively  in  a  clockwise  sense  and  ^  =  0°  along  the  minm 
axis  (b),  then  (x  =  R{^}  sin(^))„  where  R{^}  is  the  ellipse  radius  as  a  function  of  %.  The  parametn  R{^}  is  given  by 
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(9) 


[^b2sin2(^)  +  a^  cos^f^)  ^ 
and  the  sltqie,  dy/dx,  can  be  represented  as  a  fimction  of 


dy  _  -b^tanf^  (10) 

dx  ~ 


Therefore, 


a  =  90®  +  ^-tan-* 

To  maintain  a  constant  mass  flux  within  the  ellipse, 

R{’=,}  v{§}  sin(a)  =  R{§}  V,  =  CW  (12) 

where  V{^}  is  the  velocity  along  the  ellipse  perimeter  at  R{^)  as  a  fimction  of ,  Vr{^}  is  the  rotational  velocity,  and 

Cmf  is  the  mass  flux  constant  defined  by(Ciiif  =  bVb  =  aVa).  Equation  (12)  describes  the  variation  in  radius  and  velocity 
along  the  ring  of  maximum  velocity  that  separates  the  inner  and  outer  flow  regions.  Within  the  core  region,  the  rota¬ 
tional  velocity  varies  by 

vr  /  r  =  V{^}  sin(a)  /  R{^)  =  Vr{^)  /  R{^}  =  Ci{^}.  (13) 

In  the  potential  flow  region,  the  rotational  velocity  varies  by 

vr  r  =  V{^}  sin(a)  R{^}  =  Vr{^}  R{§}  =  C2{%},  (14) 

These  relationships  are  analogous  to  Eqs.  (1)  and  (2)  for  the  circular  mesocyclone. 

Examples  of  the  elliptical  mesocyclone  model  for  purely  cyclonic  rotation,  whoe  the  total  velocity  is  directed 
along  the  ellipse,  are  given  in  Figure  7.  Simulated  is  the  velocity  distribution  that  would  be  detected  by  single  Dof^ler 
radar,  with  the  radar  located  towards  the  bottom  of  the  page.  These  examples  are  for  an  arbitrary  eccentricity  of  0.866, 
which  corresponds  to  a  stretch  factor  (S)  of  2.  Stretch  factor  is  defined  as  the  ratio  of  the  ellipse  major  and  minor  axes 
(a/b).  This  means  that  that  S  =  1  for  a  circle  and  S  >  1  for  an  ellipse. 


dx 


(11) 


orientation  (degrees)  =  60 


Hgure  7.  Examples  of  The  Elliptical  Mesocyclone  Model  the  Condition  of  Pure  Rotation  (tiiat  is  divergence  is 
everywhere  equal  to  zoo).  The  velocity  distribution  that  would  be  detected  by  a  single  D(q)pla’  radar,  located 
towards  tiie  bottom  of  the  page,  is  shown.  Contour  levels  are  as  indicated  in  Figure  1.  White  circles  indicate 
the  (necise  locatitm  of  velocity  maxima.  Example  shown  is  for  an  ellipse  of  axial  stretch  factor  (S)  =  2.  El¬ 
lipse  size,  shape,  and  orientation  relative  to  the  radar  are  indicated  by  die  thick  solid  line.  The  strai^t  line  bi¬ 
secting  die  eUpse  is  oriented  along  die  major  axis.  The  exanqiles,  a,b,c,  and  d  are  for  die  same  feature 
observed  at  a  feature  orientation  angle  (9f)  of 90*,  60*,  30*.  and  0*,  respectively. 
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Figure  7.  (Continued) 


The  examples  shown  in  the  figure  represent  the  same  diecHetical  elliptical  mesocyclone  at  various  orwntatioas  to 
the  radar.  The  feature  orientation  angle  (  6f)  refers  to  the  orientation  of  the  ellipse  relative  to  tiw  radar.  6f  equal  to  0° 
means  the  minor  axis  is  aligned  with  the  radar  beam,  as  in  Figure  7d.  Positive  (negative)  values  of  ^indicate 
counterclockwise  (clockwise)  rotation  of  a  feature. 

The  figure  indicates  that  for  a  pure  rotation  scenario,  the  velocity  peaks  are  always  in  a  line  normal  to  the  radar 
beam  (^  =  0°),  regardless  of  9f.  However,  the  magnitude  of  die  detected  velocity  peaks  and  the  distance  between  them 
varies  greatly  as  a  function  of  the  ellipse  orientation  6f.  Figure  8  shows  the  effect  on  the  ai^arent  mesocyclone  shear  of 
Of  and  stretch  factor  S.  Here  the  observed  shear  is  normalized  to  the  peak  value,  which  cmesponds  toOf = 90  °. 

The  potential  impact  of  ellipse  «ientation  to  the  radar  increases  with  the  stretch  factor.  With  even  a  modest  stretch 
factor  of  1.25,  the  detected  shear  by  single-Doppler  radar  will  be  only  64  percent  of  the  true  peak  value  when  the  radar 
is  loddng  broadside  at  die  ellipse  ( Of  =  0°).  For  the  examples  of  Figure  7,  widi  S  =  2,  the  detected  shear  may  be  only  25 
percoit  of  the  peak.  This  has  inqxntant  implications  for  mesocyclone  detection.  The  established  crit«ia  for 
mesocyclone  identificatioo  by  single-Dof^ler  ra^  require  die  detection  of  ''mesocyclonic  shear''  at  several  contiguous 
elevations^^  (JDOP  Staff,  1979).  In  practice,  mesocyclonic  shear  is  defined  as  0.(X)5/s  or  greato-  across  the  relative 
inbound  and  outbound  velocity  peaks  associated  with  the  observed  velocity  couplet  Therefore,  with  die  present 
mesocyclone  detection  criteria,  the  early  detection  of  all  mesocyclones  and  the  detection  of  marginal  mesocyclones  with 
peak  shear  near  0.(X)5/s  depends  on  die  feature  orientation  relative  to  the  radar.  If  mesocyclone  elongation  and 
orientation  are  a  function  of  the  low  level  inflow,  there  will  be  prefoential  regions  for  mesocyclone  detection  within  the 
radar  sfjmning  circle.  If  the  Del  City  storm  is  representative  of  Oklahoma,  then  die  locations  fw  optimal  shear  detection 
should  be  NNW  and  SSE  of  the  radar. 

Figure  5b  of  the  Del  City  examples  is  rqnesented  furiy  well  by  the  model.  The  Del  City  example  has  a  stretch 
facta  (S)  of  about  1.9  and  a  feature  orientation  angle  (Of)  of  api»oximately  55°.  A  single-Dqppler  model  simulatioi 
with  diese  characteristics  is  displayed  in  Hgure  9  along  widi  the  Del  City  observation.  There  is  good  general  agreement 
between  the  two.  Orientation  of  dw  velocity  peaks  with  respect  to  the  radar  is  very  similar,  and  the  “zero  line”,  indicated 
by  the  white  contour  interval  (non-shaded),  has  about  the  same  orientation  with  respect  to  the  peaks,  aldiough  with  E>el 
City  the  contours  between  the  velocity  peaks  are  not  entirely  as  straight  as  in  the  model.  In  boA,  it  is  easy  to  determine 
the  ellipse  outline  dirougb  the  contour  inflections  at  the  boder  between  the  inner  core  flow  and  the  potential  flow 
region. 

The  Dei  City  example  is  at  a  height  of  1.1  km,  which  is  generally  low  for  pure  rotation  in  a  mesocyclone.  However, 
this  height  is  at  the  lower  bounds  of  where  Brandes^  found  non-divergent  flows  to  exist  for  this  particular  mesocyclone. 


3,2  Kinmnatics  of  the  Model  Rotational  Flow 

Vorticity  and  divergence  fields  of  the  model  are  determined  fa  pure  rotational  flow.  The  parameters  are  determined 
using  a  simple  point  diffoencing  formulation  from  a  Cartesian  grid.  For  pure  rotation,  divergence  is,  by  definition, 
equal  to  zero  in  botii  the  inner  and  outer  flow  regions  of  the  ellipse.  The  voticity  field  is  more  interesting  and  is  shown 
in  Figure  10.  As  with  the  circular  flow,  vorticity  is  uniform  wiAin  the  ellipse  coe  region  and  is  equal  to  the  average  of 
the  radial  shear  along  the  maja  and  mina  axes. 

The  figure  shows  examples  where  the  inner  vorticity  is  held  constant  while  the  stretch  facta  (S)  is  varied  between 
3.0, 2.0, 1.5,  and  1.0.  In  the  outer  region,  the  elliptical  flow  produces  a  pair  of  anticyclonic  vorticity  peaks  located  either 
side  of  the  ellipse  mina  axis.  The  magnitude  of  this  negative  vorticity  relative  to  the  positive  vorticity  of  the  cote  region 
varies  as  a  fiin^on  of  die  eccentricity;  the  greater  the  stretch  factor,  the  greater  die  relative  anticyclonic  component 

While  dw  couplet  of  anticyclonic  vorticity  is  a  by-product  of  the  model,  does  it  have  some  validity  in  real  wald? 
The  examples  in  Figure  7  are  fa  a  rotational  wind  field  that  varies  according  to  a  general  potential  flow  in  die 
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Normalized  Shear 


Observed  shear  as  function  of  ellipse  orientation 


ellipse  orientation  (degrees) 


Hgure  8.  Detectable  Shear,  by  Single-Doppler  Radar,  as  a  Function  of  the  Mesocyclone  (ellipse)  Orientation  with  Re¬ 
spect  to  the  Radar.  The  effea  of  stretch  factor  is  examined  with  S  varying  from  1.0  (circular  flow)  to  3.0. 
Peak  shear  is  detected  when  6f  =  90°  for  pure  rotation  (no  divogence).  Other  shears  within  each  category  of 
S  are  normalized  to  this  value. 
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DEL  CITY 


Comparison  of  the  Model  Output  for  Pure  Rotation  (shown  in  a)  to  a  Real  Life  Example  From  the  Del  City  Storm  (shown  in  b).  Hie  model  example 
approximates  the  Del  City  mesocylone  size  and  shape  characteristics,  and  orientation  to  the  radar.  Presented  scales  of  size  are  different.  Contour  lev¬ 
els  are  as  indicated  in  Figure  1. 


r 


axial  stretch  factor  =3.0 


Hgure  10.  Presentations  of  the  Vcrticity  Held  for  an  Elliptical  Mesocyclone  Flow  as  a  Function  of  the  Stretch  Factor. 

Peak  velocities  vary  for  each  example  to  maintain  the  same  vorticity  magnitude  in  the  core,  0.02S/s.  The 
feature  is  presented  at  a  orientation  angle  of  9f = 4S’,  to  show  die  vorticity  contours  as  clearly  as  possible, 
but  the  vorticity  field  is,  of  course,  not  a  function  of  the  feature  orientation  to  the  radar.  The  exanqile  show 
how  the  peak  in  andcyclonic  vcrticity  outside  the  core  varies  as  a  function  of  stretch  factor;  a)  S  =  3.0,  b) 
S  =  2.0,c)S=  1.5,d)S=  1.0. 
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d) 


Hguie  10.  (Contined) 
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aiesocyclone  outer  flow  region,  as  described  in  Eq.  (14).  While  true  potential  flow  is  unlikely,  it  is  intnesting  to 
examine  how  the  flow  in  the  outer  region  might  be  varied  to  eliminate  the  anticyclonic  couplet  entirely.  Figure  1 1  shows 
how  the  velocity  field,  seen  by  single-Doppler  radar,  would  have  to  vary  to  have  no  associated  anticyclonic  vorticity 
couplet.  Such  elongation  of  the  velocity  contours  in  the  outer  flow  region  is  not  seen  in  single-Dof^ler  mesocyclone 
observations,  and,  therefore,  the  presence  of  negative  vorticity  pairs  in  a  general  elliptical  mesocyclone  flow  is  at  least 
dynamically  possible.  This  is  interesting  considering  that  some  cyclonic  tornadoes  form  in  conjunction  with  an 
anticyclonic  counterpart  that  is  generally  of  lesser  intensity.  Whether  such  an  anticyclonic  vorticity  zone  serves  as  the 
source  for  amplification  to  tornado  intensity  is  an  important  question,  but  beyond  the  present  scope  of  this  work. 


33  Rotattonal  Kinetic  Energy 


A  convenient  parameter  for  describing  mesocyclone  organization  is  Rotational  Kinetic  Energy  RKE.  While 
mesocyclone  RKE  is  not  especially  effective  for  evaluating  tornado  potential^^  (Donaldson  and  Desrochers,  1990),  it  is 
a  conservative  quantity  that  may  be  useful  for  evaluating  other  types  of  severe  weather,  and  is  worthy  of  discussion 
here. 


It  is  instructive  to  derive  the  equation  for  circular  RKE.  For  circular  mesocyclones  it  is  convenient,  and  reasonable, 
to  assume  solid  body  rotation  where  vorticity  is  uniform  in  the  core.  Rotational  kinetic  energy  is  then  simply  (I  aP-/2), 
the  product  of  moment  of  inertia  (I)  and  angular  velocity  (a>).  For  a  solid  disk,  (I  =  m  R^/2),  the  product  of  the 
mesocyclone  mass  (m)  and  radius  (R)  of  the  core.  Mass  is  given  by  (p  Ah  n  R^).  Angular  speed  (to)  is  given  by  Vi/R  . 
Rotational  kinetic  energy  for  a  circular  mesocyclone  is  represented  by: 

RKEcirele  =mV,^/4  =  pAhJcR2V,^/4  (15) 


The  identical  solution  can  be  obtained  by  integrating  the  kinetic  energy  of  a  sector  around  the  entire  mesocyclone 
circle.  It  is  necessary  to  take  this  jqiproach  for  an  elliptical  mesocyclone  because  the  flow  cannot  be  described  as  solid 
body,  although  the  vorticity  field  inside  the  core  may  be  uniform.  Also,  an  elliptical  flow  does  not  have  a  moment  of 
inertia  since  the  elliptical  shrqre  is  not  rotating  around  itself,  but  rather  the  flow  is  following  a  relatively  steady  elliptical 
sh^.  We  begin  with  the  rotational  kinetic  enogy  of  an  ellipse  sector,  q  radians  wide, 

RKE  s  =  ms  (Vr  ©cm  {^})^/2  (16) 

where  Vr  ecm{^}  is  the  rotational  velocity  at  the  sector  center  of  mass  as  a  function  of  the  ellipse  relative  angle  The 
mass  of  a  sector  is  given  by  ms  =  p  Ah  A$.  The  area  of  a  sector  is  given  by  (As  =  R{^)^  dri/2  =  Rem  (4}^  dq),  where  dq 
is  the  sector  angular  width,  assumed  small,  and  Rcm{4}is  the  radius  to  the  sector  center  of  mass.  The  velocity  along  the 
ellipse  at  Rcm{^}  is  Vcm{^}  and  is  given  by(Vcm{^)  =  o>{^}  Rcm{^}  =  V{^}  Rcm{^}  /R{^})  •  Sirjce  the  rotational  ve¬ 
locity  is  a  function  of  a;  (Vr  9cm{^}  =  Vcmli)  sin(o)  =V{^)  sin(a)  Rcm{^}  /R{^}).  From  Eq.  (12),  (Vj^}  sin(a)  =  aV 
a  /  R{^}),  and  substituting  for  Rem  {^}>  the  rotational  kinetic  energy  fOT  an  elliptic  sector  is  given  by 

RKEs  =  pAha2(Va)2dq/8.  (17) 

Integration  of  q  from  zero  to  2  n  to  obtain  the  total  ellipse  rotational  kinetic  energy  yields: 

RKEeUipse  =pAha2(Va)27t/4  =  mea(Va)2/4b  (18) 

where  the  ellipse  mass  (me)  =  p  Ah  n  a  b.  For  a  circle,  where  a  =  b,  Eq.  (18)  reduces  to  Eq.  (IS). 
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rotatlcMi  (degrees)  =  45 
axial  stretch  factor  s  2. 0 


rotation  (degrees)  =  90 
axial  stretch  factor  =2.0 


Figure  11.  An  Example  of  How  the  Velocity  Field  in  the  Potential  Flow  Area  of  an  Elliptically-shaped  Mesocyclone 
Would  Have  to  Vary  to  Have  Zero  Anticyclonic  Vorticity.  Contour  levels  are  as  indicated  in  Figure  1 . 
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ellipse.  Therefore,  ellipse  RKE  can  be  determined  without  knowledge  of  the  eccentricity,  ellipse  orientatitm  to  the  radar, 
or  the  magnitude  of  Va  or  Vb.  The  single-Doppler  velocity  peaks  can  be  sf^lied  to  the  relationship  for  a  circular 
mesocyclone  Eq.  (IS)  to  yield  the  RKE  of  the  ellipse,  regardless  of  the  eccentricity.  Mesocyclone  rotational  kinetic 
energy  is  dmefore  a  conservative  quantity  for  detection.  Unfortunately,  in  our  own  investigations,  RKE  has  not  been 
found  to  be  of  significant  value  for  evaluating  the  mesocyclone  for  hail  or  tornado  potential.*^’*^ 


3.4  Excess  Rotational  Kinetic  Energy 


Donaldson  and  Desrochers’^  present  the  concept  of  Excess  Rotational  Kinetic  Energy  (ERKE),  which  has 
subsequently  been  shown  to  offer  [xoinise  for  hail  and  especiaUy  tornado  prediction  (for  example,  Desrocbos  and 
Donaldson,^).  From  their  work,  ERKE  for  a  circle  is  given  by 


ERKEcircle  =  pAhnR2(Vr-RSm)2/4  (19) 

where  Sm  is  an  arbitrary  shear  threshold  that  defines  mesocyclone  rotation,  nominally  O.OOS/s  (after  Donaldson^).  ERKE 
is  defined  to  exist  only  for  (Vr  -  R  Sm)  >  0.  ERKE  is  simply  the  rotational  kinetic  energy  in  excess  of  what  would  be 
needed  to  sustain  mesocyclone  rotation.  The  idea  of  ERiS  is  to  distinguish  large,  well-organized,  and  swiftly  rotating 
features  that  may  present  a  severe  weather  threat  from  feamres  that  may  have  greater  RKE  but  are  benign. 


ERKE  fn-  an  ellipse  can  be  determined  in  an  analogous  manner  to  RKEeiiipce.  The  ERKE  of  an  ellipse  sector  is 
represented  by 

ERKEs  =  p  R{4}2  (V{^}  sm(a)  -  R{^}  Sm)^  d  r\/  8.  (20) 


Here,  in  analogy  with  the  circular  flow,  ERKE  exists  for  (V{^}  sin(a)  >  R{^}Sni)-  Since  the  radial  shear  varies  accord¬ 
ing  to  the  ellipse  relative  angle  (^.  it  is  possible  that  radial  shear  satisfies  Sm  for  only  a  sector  of  the  ellipse.  Substituting 
Eq.  (9)  for  R{^}  into  the  inequality  and  solving  for  the  maximum  ellipse  relative  angle  fix  integrating  Eq.  (20)  for 
ERKE  in  the  fint  quadrant  of  the  ellipse  is 


Expansion  and  integration  of  Eq.  (20)  for  each  ellipse  quadrant  yields  the  ERKE  for  an  ellipse: 


(21) 


ERKEempse  =  B^  /  (^a2 Va2-2a V,SmR(4f +  R{^rSm2)di 
4  =  0^ 


(22) 


Figure  12  shows  Eq.  (22)  evaluated  using  Simpson’s  rule  and  normalized  by  the  corresponding  value  of  RKE. 
Here,  the  stretch  factor  (S)  is  varied  between  1.0  and  3.0,  and  die  shear  threshold  Sm  is  varied  between  0.001/s  and 
0.009/s.  RKEeiiipse  is  kept  constant  for  the  calculations  and  is  determined  for  a  climatological  mature  mesocyclone  of 
2.75  km  radius  and  23.3  m/s  rotational  velocity  (after  Burgess  et  al.’^). 

With  RKEeiiipse  kept  constant  as  a  reference,  ERKE  decreases  as  the  mesocyclone  becomes  more  elongated, 
because  the  average  shear  decreases.  The  change  with  Sm  is  a  function  of  the  products  (a  Va)  and  (b  Vb),  which  are 
proportional  to  RKE.  For  a  low  shear  threshold  (Sm  =  0.001/s),  ERKE  decreases  by  26  percent  between  S=1  and  S=3. 
At  progressively  higher  shear  thresholds,  the  relative  reduction  in  ERKE  increases.  For  example,  at  Sm  =  O.OOS/s  ERKE 
is  diminished  by  about  72  percent  between  S=1  and  S=3. 
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Figure  12.  Variation  of  ERKE  as  a  Function  of  Ellipse  Stretch  Factor  (S)  and  ERKE  Shear  Threshold  (Sm).  Values  are 
normalized  by  RKE  calculated  for  b  =  2750  m  and  Vb  =  23.3  m/s,  climatological  values  of  an  average  ma¬ 
ture  mesocyclone  (from  Burgess  et  al.^^). 


Another  way  to  examine  ERKE  is  to  cooq>aie  it  to  vorticity.  Figure  13  shows  ERKE  vahatitHi  widi  vorticity  and 
stretch  hKtor.  Here,  ERKE  at  Sm  =  0.005/s  is  calculated  for  various  vorticity  levels  and  is  normalized  to  the  ERKE 
value  at  S  -  3.  For  any  given  magnitude  of  mesocyclonic  vorticity  widiin  die  core  region,  ERKE  increases  as  a  feature 
becomes  mme  elongated.  The  relative  change  of  ERKE  decreases  with  increasing  vwticily. 

While  the  increase  in  ERKE  widi  feature  elongation  at  constant  vorticity  is  notable,  it  would  not  be  necessarily 
correct  to  infer  that  ERKE  would  indicate  a  significant  alteration  in  analysis,  that  is  in  storm  strength  «  structure,  based 
solely  on  elongation.  In  fact,  the  variation  in  ERKE  widi  elongadcm  is  small  conqiared  to  the  variation  urith  vorticity. 
The  Del  City  example  of  Figure  Sb  is  an  aiqiropriate  case  to  make  this  point  The  time  of  this  particular  scan  was  only  S 
min,  before  the  first  F2  tornado.  Vrvticity  here  is  estimated  at  0.02/s,  and  with  S  =  1.9,  the  mesocyclrme  was 
significandy  elongated.  Crunparing  the  ERKE  of  this  feature  to  die  ERKE  of  a  rntne  elongated  feature  at  a  lesser 
vorticity,  and  to  a  less  elongated  feature  at  a  higher  voticity  gives  an  indication  of  the  importance  of  vorticity  in  ERKE. 
ERKE  of  the  Del  City  case  is  roughly  17  times  greats  in  magnitude  than  ERKE  fm  a  feature  of  one-half  die  vorticity 
(0.01/s)  and  S  =  3.  Experience  indicates  that  such  a  maigiiial  vorticity  case  would  ncrt  be  eiqiected  to  produce 
significant  tornadoes  (for  example,  see  Desrochers  and  Donaldson^).  Con^iare  this  to  a  mesocyclone  of  twice  die 
vorticity  of  die  Del  City  mesocyclone.  ERKE  for  a  feature  of  S  =  1,  a  circular  case,  and  vorticity  0.04/s,  is  roughly  4 
times  that  of  Del  City.  This  is  in  the  range  of  what  is  observed  for  violent  tornado  producing  strains.  For  each  doubling 
in  vorticity  magnitude  there  is  nearly  an  order  of  magnitude  increase  in  ERKE,  which,  excqpt  for  very  marginal 
mesocyclones,  is  greatra  than  die  change  attributable  to  eccentricity.  This  assumes  that  S  =  3  is  a  safe  upper  bound  fra 
eccentricity.  In  summary,  vorticity  contributes  more  to  die  variation  in  ERKE  than  eccentricity.  At  die  marginal 
mesocyclone  vorticity  level,  eccentricity  is  perhaps  of  litde  significance  because  it  can  not  contribute  to  significant 
values  of  ERKE.  At  increasingly  larger  vraticity  values,  where  die  absolute  magnitudes  of  ERKE  are  already  significant, 
it  is  speculated  diat  increased  eccentricity  may  add  to  the  apparent  overall  severity  of  die  mesocyclone. 

The  investigation  of  ERKE  in  mesocyclones  has  to  date  assumed  that  the  mesocyclone  flow  is  circular  (for 
example,  Donaldson  and  Desrochers^^).  Given  the  considerable  variability  in  the  apparent  size  and  rotational  velocity  of 
an  elliptical  mesocyclone,  it  is  interesting  to  ask  fra  what  feature  orientation  a  “circular”  assumption  would  be  accurate. 
In  other  words,  at  what  orientation  of  the  ellipse,  relative  to  die  radar,  will  the  observed  velocity  peaks  be  of  sudi 
magnitude  and  azimuthal  separation  that  the  Excess  Rotational  Kinetic  Energy  (ERKE)  of  the  ellipse  can  be  draived 
under  the  assumption  that  the  peaks  are  associated  with  a  circular  flow? 

Figure  14  shows  the  variation  of  “apparent”  ERKE  for  an  elliptical  flow  where  the  vorticity  of  the  core  and  die 
shear  threshold  (Sm)  are  kept  constant  ^pses  of  various  stretch  factors  are  examined  and  viewed  from  diffnent 
orientations.  The  vorticity  magnitude  O.OOT/s  has  been  selected  to  be  that  for  climatological  mature  mesocyclones.  The 
shear  threshold  is  arbitrary  set  at  0.(X)S/s.  The  velocity  field  is  rate  of  pure  rotation.  ERKE  is  determined  assuming  the 
radar  observed  velocity  peaks  are  associated  with  a  circular  feature.  These  values  are  dien  normalized  by  the  true  ERKE 
of  the  feature. 

While  the  results  correspond  to  a  climatological  feature,  die  generalities  revealed  are  important  ERKE  estimation 
errors  will  be  minimized  when  a  feature  is  oriented  about  ±  45°  to  the  radar  beam.  For  orientations  at  larger  angles, 
ERKE  may  be  overestimated  by  a  factor  of  2  or  mrae  as  the  maximum  velocities  along  the  ellipse  minor  axes  are 
detected.  At  orientations  less  than  45°  ERKE  will  be  underestimated,  so  sevraely  that  a  severe  mesocyclone  may  tqipear 
to  have  zero  ERKE.  The  amount  of  over-  or  undraestimation  will  increase  as  a  feature  becomes  more  elongated.  Even 
with  modest  eccentricity,  the  errors  can  be  considerable.  With  a  stretch  factra  of  S  =  1.2,  ERKE  estimation  can  be  off  by 
as  much  as  SO  percent  if  the  eccentricity  is  not  known. 

4.  DIVERGENCE  SCENARIOS  FOR  ELUPTICAL  FLOW 


The  divergence  field  for  a  circular  flow  has  been  discussed  in  Section  2.  In  this  section,  a  few  possible  divergence 
scenarios  fra  an  elliptical  flow  are  examined. 
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Figure  13.  Variation  of  ERKE  with  Feature  Eccentricity  and  Vorticity.  Values  of  ERKE  for  each  given  vorticity  level 
are  nnmalized  to  the  ERKE  at  a  stretch  factor  of  3.0. 
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Figure  14.  Hie  Ratio  of  Circle  ERKE  to  the  True  Ellipse  ERKE  as  a  Function  of  Stretch  Factor  (S)  and  Feature  Orien¬ 
tation  (0f).  The  intent  is  to  illustrate  potential  errors  in  ERKE  estimation  by  qiplication  of  the  old,  circular 
mesocyclone  model.  Circle  ERKE  is  determined  from  the  single-Doppler  mesocyclone  observation,  inter¬ 
preted  to  be  associated  with  a  circular  flow,  regardless  of  how  elliptical  the  true  flow  actually  is.  Vorticity  is 
held  constant  for  these  calculations. 


25 


4.1  Untform  DhFWfaice 


A  logical  starting  point  is  to  sinq>ly  reapply  die  mesocyclone  model  for  divergence,  such  that  flow  is  directed 
normal  to,  instead  of  idtmg,  the  ellipse.  Namely,  (R{^}  V(4}  sin(a)  =  R{^}  V(i{^}  =  a  Va),  where  Vd{4}  is  the 
divergent  component  of  the  velocity.  A  change  in  the  dir^tm  of  flow  m  this  manner  results  in  a  flow  field  that  has  zero 
vofticity  everywhere,  and  uniform  divergence  in  the  core  region.  Mesocyclone-like  signatures  in  the  Dof^ler  field  are 
not  observed  for  purely  convergent  signatures,  but  signatures  indicative  of  a  purely  divmgent  flow  are  sometimes  seen  at 
die  iqiper  teaches  (tf  a  storm.  A  model  example  of  the  divergent  field  is  shown  in  Figure  IS,  which  is  analogous  to  the 
model  vordcity  field  (Hgure  10b).  Uniform  divergence  exists  in  the  mesocyclone  con  region  and  extends  into  die 
potential  flow  region  along  the  major  axis.  Along  the  minor  axis,  in  the  potential  flow  rpgicm,  a  divergence  couplet  of 
opposite  sign  to  that  of  the  core  r^on  exists. 

An  exariqile  of  the  single-Doppler  rqipearance  of  the  velocity  field  associated  with  the  divergence  field  of  Figure 
IS  is  shown  in  Figure  16  at  various  radar-relative  ellipse  orientation  angles.  These  examples  are  meant  to  represent 
upper  storm  outflow,  as  can  be  noted  by  the  signs  of  the  velocity  contours.  The  flow  is  characterized  by  opposing 
velocity  peaks  oriented  in  the  radial  direction,  =  90°).  Straight  isodops  in  the  core  region  mean  uniform  divergence 
exists  thm.  Although  die  magnitudes  of  the  peak  velocities  change  as  the  ellipse  is  rotated  with  respect  to  the  radar,  the 
orientation  of  the  peaks  remains  die  same,  as  was  also  observed  for  pure  rotation  (Figure  8).  It  is  interesting  to  see  in 
Figure  15,  that  in  a  general  divergent  flow  thoe  is  a  convergence  couplet  in  the  potential  flow  region.  This  couplet  is 
located  along  the  direction  of  the  minor  axis  where  the  outflow  velocity  is  maTimum  This  result  is  produced  by  the 
assunqition  that  the  velocity  magnitnA*  stricdy  adheres  to  a  potential  flow  rate  of  change,  as  in  Eq.  (14),  which  may  not 
be  realistic. 


4.2  Urdfonn  Convergeiice  ruid  Rotation 


At  most  elevations  through  a  mesocyclone,  the  flow  field  is  conqirised  of  a  combination  of  divergence  and 
vOTticity.  Hoe  we  examine  low  level  inflow.  Figure  17  presents  single-Doppler  views  when  die  magnitude  of 
convergence  in  the  core  region  is  equal  to  diat  of  vorticity.  This  is  shown  for  an  arbitrary  stretch  factor  (S)  of  2. 

The  white  dots  in  the  figure  indicate  the  location  of  die  maximum  radar-relative  velocities.  The  orientation  of  die 
velocity  peaks  relative  to  the  radar  does  not  change  if  the  feature  is  observed  fiom  different  angles.  Therefore,  when 
divergence  and  vorticity  are  uniform  within  die  core,  as  in  dwse  examples,  the  ratio  of  divergence  to  vorticity  is  a 
function  of  the  velocity  peaks’  mientation  angle  (4): 


Divergence 

Vorticity 


=  tan  (<►), 


(23) 


as  is  also  true  for  the  circular  model  (Section  2).  This  is  a  very  useful  quality,  although  die  eccentricity  must  be  known 
to  determine  the  absolute  magnitudes  of  divergence  and  vorticity. 

Figure  Sa,  presented  again  for  convenience  in  Hgure  18b,  shows  an  example  fiom  Del  City  where  the  elliptical 
mesocyclone  core  rqipears  to  be  comprised  of  a  combination  of  rotation  and  convergence.  A  velocity  peak  (mentation 
angle  ^  of  31°  indicates  a  divergence/vorticity  ratio  of  -0.6.  This  relative  magnitude  of  convmgence  in  the  flow  is 
consistmit  with  a  0.5  km  height 


The  model  simulaticm  (Figure  18a)  compares  well  with  the  real-life  example.  The  “zero  line”  of  both,  and  the 
orientation  of  the  vel(x:ity  peaks,  are  similar.  The  relative  straightness  of  the  velocity  contours  widiin  die  core  region  of 
the  Dei  City  exanqile  would  suggest  that  both  convergence  and  vorticity  are  nearly  uniform  within  the  core.  Velocity 
contours  of  the  “potential  flow”  region  do  not  agree  as  well.  For  one  thing,  there  appears  to  be  a  change  in  the  potential 
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axial  stretch  factor  =2.0 


Hgure  IS.  Exaiiq>le  of  the  Divogence  Field  Produced  by  the  EUiptical  Mesocyclone  Model.  Divergence  is  uniform 
within  the  core.  The  divotgence  field  is  indq)endent  of  the  fixture  orientation  presented.  Divergence  for  a 
stretch  factor  of  2  is  shown. 


Hgine  16.  Examples  of  the  Elliptical  Mesocyclone  Model  for  Pure  Divergence,  y^diin  the  ctne,  the  velocity  field  cor¬ 
responds  to  uniform  divergence.  Outside  the  cote,  velocity  tapos  off  according  to  a  standard  potential  flow. 
The  viewpoint  is  diat  of  a  single-DopplCT  radar,  located  towards  the  bottom  of  the  page.  Contour  levels  are 
as  indicate  in  Hgure  1.  White  circles  indicate  the  precise  locations  of  velocity  maxima.  Example  shown  is 
for  an  ellipse  of  stretch  factor  (S)  =  2.  Ellipse  size,  shape,  and  orientation  relative  to  die  radar  are  indicated 
by  die  thick  solid  line.  The  straight  line  bisecting  die  ellipse  is  along  the  majtv  axis.  The  examples,  a,b,c, 
and  d  are  fOT  the  same  feature  obsoved  at  a  feature  cnientadon  angle  (9  f )  of  90,  60,  30,  and  0*,  respec¬ 
tively. 
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Figure  17.  E3rainples  of  the  Elliptical  Mesocyclone  Model  Velocity  Field  for  Equal  Magnitudes  of  Convogence  and 
VoTticity  Inside  the  Core.  The  viewpoint  is  that  of  a  singJe-Doppler  rato,  located  towards  the  bottom  of  die 
page.  Contour  levels  are  as  indicated  in  Figure  I.  White  circles  indicate  the  fnecise  location  of  velocity 
maxima.  Example  shown  is  for  an  ellipse  of  stretch  factor  (S)  =  2.  Ellipse  size,  shape,  and  mentation  rela¬ 
tive  to  the  radar  are  indicated  by  the  thick  solid  line.  The  straight  line  bisecting  the  ellipse  is  along  the  ma¬ 
jor  axis.  The  examples,  a,b,c,  and  d  are  for  the  same  feature  observed  at  a  feature  orientation  angle  (9f )  of 
90, 60, 30,  and  O",  respectively. 
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Figure  18.  Comparison  of  the  Elliptical  Mesocyclone  Model  a)  to  an  example  from  the  Del  City  storm  (b).  Each  show  a  flow  field  comprised  of  a  combination 
of  convergence  and  rotation.  The  Del  City  example  is  estimated  to  have  a  stretch  factor  (S)  of  1.9,  is  oriented  at  6f  =:  SS°  and  convergence  within 
the  core  is  equal  to  0.6  of  the  vorticity.  The  model  example  represents  these  same  characteristics.  Presented  scales  of  size  are  somewhat  different. 
Contour  levels  are  as  indicated  in  Figure  1 . 


flow  from  one  side  of  the  mesocyclone  to  the  other,  which  is  not  included  in  the  model.  Velocities  on  the  right  side 
(figure  relative)  of  the  Del  City  mesocyclone  change  at  a  rate  greater  than  potential  flow.  The  0{qx>site  is  true  of  the  left 
side.  Clearly,  in  the  exaiiq>les  seen  up  to  now,  only  the  core  region  fits  our  simple  model. 


4 J  Non-Uniform  Convergmce 


In  previous  sections,  the  investigation  has  addressed  conditions  where  uniform  vorticity  anj  divergence  fields 
exist  in  the  mesocyclone  core  region.  In  this  section,  the  condition  of  non-uniform  fields  is  addressed.  Non-uniform 
convergence  is  likely  to  exist,  fm  example,  at  low  levels  of  a  stcvm  where  the  inflow  is  unidirectional.  Natural  variations 
in  inflow  patterns  are  likely  irmumerable.  It  is  helpfiil  to  examine  cases  where  the  vorticity  and  divergence  fields  are 
known. 

In  the  Brandes  (1984)  multiple-Doppler  study,^  aspects  of  the  flow  of  the  Del  City  storm  are  revealed.  Figure  19a 
shows  his  analysis  of  the  convergence  term  contributing  to  vorticity  tendency  at  1.3  km  AGL  and  1826  CST  (here,  north 
is  towards  the  top  of  the  figure).  Also  shown  is  the  vertical  vorticity  at  the  same  height  (Figure  19b).  The  dashed  line  in 
each  figure  indicates  the  mesocyclone  outline  as  defined  by  vorticity  >  0.01/s.  The  thick  solid  line  is  the  mesocyclone 
outline  determined  through  single-Dc^pler  analysis.  Both  of  these  indicate  an  elliptical  flow  pattern.  In  this  example  of 
the  Del  City  mesocyclone,  vorticity  was  non-uniform  in  the  “core”  region,  peaking  inside  the  cote  in  a  region  southeast 
of  the  geometric  center. 

A  perspective  of  the  inflow  is  obtained  from  the  convergence  term  (Figure  19a),  given  by  (  dux/dx  +  3uy/9y) 
(duy/dx  -  dux/dy).  The  convergence  term  of  the  vorticity  tendency  equation  is  not  necessarily  representative  of  the  wind 
field  convergence.  However,  if  vorticity  is  uniform  the  convergence  term  will  indicate  the  distribution  of  divergence, 
although  not  the  magnitude.  Figure  19b  shows  that  vorticity  is  not  uniform  in  the  core.  However,  this  may  be  die  result 
of  a  non-uniform  inflow,  since  the  vorticity  peak  is  off-center  from  the  mesocyclone  ellipse.  This  contention  will  be 
supported  if  the  shape  and  magnitude  of  the  convergence  term  seen  in  Figure  19  a  and  b  can  be  simulated  by  including 
rotation  derived  from  the  elliptical  model  that  is  centered  at  the  ellipse  geometric  center  and  contains  a  constant 
horizontal  mass  flux. 

The  convergence  field  was  simulated  by  trial  and  oror,  with  guidance  from  the  horizontal  wind  vectors  (Figure 
19c),  which  indicate  convergence  on  the  mesocyclone’s  east  side.  The  most  satisfactory  reproduction  is  achieved  when 
die  maximum  inflow  is  set  at  ^  =  50°  and  the  minimum  set  at  ^  =  160°,  with  the  mesocyclone  minor  axis  (^  =  0°) 
oriented  towards  the  NE.  Furtheimcxe,  at  ^  =  160°  the  wind  field  is  slighdy  divmgent  such  that  the  peak  velocity  due  to 
divergence  is  10  percent  of  the  peak  inflow  value  at  ^  =  50°.  Therefore,  the  incoming  velocity  normal  to  the  ellipse 
surface  is  approximated  by  the  relationship: 

Vd  { 160°  <  ^  <  50°}  =  a  Va cos  (I  §  -  50°l  /1. 15)  /  R{^}.  (24) 

For  other  angles, 

Vd  {0°  <  ^  5  50°}  =  a  Va  cos  (^.61)  /  R{^},)  (25) 


Vd  { 160°  <  ^  S  50°}  =  a  Va  cos  ((^  +  50°)  /2.61)  /  R{^}.  (26) 

This  relationship  is  applied  from  the  mesocyclone  geometric  center,  and  velocity  magnitudes  vary  linearly  along  any 
particular  radial 
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Figure  19.  Fields  of  the  Convogence  Term  of  the  Vorticity  Tendency  Equation  (a).  Vertical  Vorticity  (b),  and  Low 
Level  V^ds  (c)  Determined  From  Multiple-doppler  Analysis  of  the  Del  City  Storm  at  1S26  CST  and  1.3 
km  AGL  (ad!q>ted  from  BrandesJ).  The  convergence  term  has  units  of  10^/s^.  Vorticity  units  are  10^/s.  In 
(c),  horizontal  wind  vectors  are  shown  with  radar  reflectivity  superinqwsed,  and  the  stippled  area  towards 
the  top  of  the  figure  indicates  the  region  of  subsequent  tornado  damage  that  began  14  minutes  after  this 
time.  The  medium-thickness  solid  line  reproduced  in  each  figure  is  the  30  dBz  reflectivity  contour.  The 
dashed  line  in  each  is  the  objectively  determined  mesocyclone,  corresponding  to  the  0.01/s  vorticity  con¬ 
tour.  The  thick  solid  fine  surrounding  this  is  the  mesocyclone  outline  determined  fr^om  subjective  analysis 
of  the  single-Doppler  velocity  field. 


The  simulated  convergence-teim  field  fm  the  mesocyclone  area  of  Figure  19a  is  shown  in  Figure  20.  This 
represents  an  ellipse  of  stretch  factor  (S)  =  2.0  with  the  peak  inflow  and  rotational  velocities  adjusted  to  ap(m>ximate  the 
convergence  term  field  of  the  multiple-Dq>pler  analysis.  As  in  Figure  19,  the  mentation  of  the  field  is  presented  to 
correspond  to  its  geogr^riiic  orientation  and  not  to  the  radar.  Although  the  simulated  field  is  not  as  detailed  as  the  true 
field,  it  is  representative  in  magnitude  and  general  appesnace. 

The  corresponding  divergence  field  is  shown  in  Figure  21.  Convergence  in  the  core  region  is  greatest  on  the 
mesocyclone’s  eastern  side,  and  there  is  slight  divergence  in  the  southwestern  quadrant.  One  aspect  that  is  not  well 
represented  is  the  region  of  potential  flow.  This  is  especially  ^iparent  in  the  ncxtheast  quadrant  where  the  flow  is 
divergent  No  attempt  is  made  here  to  improve  the  representation  of  die  potential  flow  since  the  immediate  results 
concern  only  the  flow  inside  the  core. 

One  effect  of  a  non-uniform  convogence  field  in  the  core  is  an  induced  vorticity  field,  shown  in  Figure  22.  The 
resultant  vc^city  field  from  the  combined  fields  of  induced  vorticity  and  vorticity  of  pure  rotation  is  shown  in  Figure 
23.  In  this  simulation,  vorticity  is  largest  in  the  southeast  quadrant  of  the  mesocyclone  core,  in  general  agreement  with 
the  muldple-Doppler  derivation  (Figure  19b).  The  simulation  underestimates  the  peak  in  vmticity  by  IS  percent, 
however. 

Also  represented  in  the  simulation  is  the  region  of  anticyclonic  vorticity  in  the  potential  flow  region  to  the 
southwest  of  the  mesocyclone  core.  Little  significance  is  attached  to  this  since  the  potential  flow  region  is  not  modeled 
well.  It  can  be  argued  that  the  vorticity  distribution  in  the  potential  flow  region  may  be  due,  in  part,  to  the  inflow  field, 
since  that  is  what  is  observed  in  the  core.  However,  at  higher  elevations  of  die  Del  City  storm  at  this  same  time,  there  is 
evidence  in  the  single-Do{^ler  observations  that  a  rear-flank  downdraft  is  forming,  characterized  by  an  anticyclonic 
velocity  couplet  (not  shown).  A  developing  rear-flank  downdraft  may  also  have  contributed  to  die  high  value  of 
anticyclonic  vorticity  seen  in  the  multiple  Doppler  analysis.  Thoe  is  no  evidence  in  the  multiple  Doppler  analysis  to 
support  the  region  of  anticyclonic  vorticity  indicated  by  the  model  to  the  northeast  of  the  mesocyclone  core. 

The  single-Doppler  B-scan  view  of  the  multiple-Doppler  analysis  is  shown  in  Figure  24  along  with  the  model 
simulation.  This  is  a  radar  perspective  from  the  west  of  the  mesocyclone.  The  thick  line  in  the  model  indicates  primary 
inflection  points  along  the  velocity  contours  and  outlines  the  mesocyclone  core,  with  a  stretch  factor  of  sqqnoxinutely 
2.0.  One  distinctive  character  of  this  exanq)le  is  the  shifting  of  the  velocity  peaks  so  that  they  lie  at  one  end  of  the  cote 
region,  a  result  of  the  increased  convergence  and  vmticity  here.  The  model  simulation  is  similar.  All  the  major 
conqranents  of  the  flow  are  represented,  fix>m  the  shifting  of  the  velocity  peak  locations  to  the  elongation  of  the  contours 
on  the  “right”  side  of  die  mesocyclone.  The  one  characteristic  of  the  simulation  that  is  not  well  represented  is  the 
bending,  or  “S”  shtqie,  of  die  zero-line  seen  in  the  real-life  example.  The  generally  straight  zero-line  of  die  simulation  is 
a  result  of  rqiplying  a  linear  variation  of  velocity  in  the  core  region  along  a  given  radial  It  is  possible  to  maintain  a 
constant  hmizontal  mass  flux  for  rotation  with  veloci^  varying  other  dian  linearly.  The  result  would  be  that  vorticity  in 
the  core  would  not  be  uniform  but  a  ring  of  maximum  vcnticity  would  exist  somewhere  within  the  ccne  or  could  peak 
near  the  mesocyclone  center.  This  scenario  would  account  for  the  S  pattern  seen  in  the  real  data  and  would  explain  why 
vorticity  was  underestimated  in  the  simulation.  It  also  explains  why  the  mesocyclone  size  by  vorticity  magnitude  from 
Brandes^  is  significantly  different  from  the  single-Doppler  estimation  of  the  core  (Figure  19b). 

In  summary,  a  directional  and  non-uniform  inflow  field  was  superimposed  upon  an  otherwise  uniform  rotation  field 
that  is  describ^  by  a  simple  elliptical  model.  The  non-uniform  inflow  produced  non-uniform  divergence  in  the 
mesocyclone  core  region,  and,  mote  to  the  point,  contributed  to  the  vorticity  field,  making  it  non-uniform  as  well. 

5.  DISCUSSION  AND  CONCLUSIONS 


To  the  author’s  knowledge,  all  single-Doppler  work  of  mesocyclones  to  date  has  been  based  on  the  assumption  that 
the  mesocyclone  can  be  iq>proximated  by  a  circular  flow.  This  wotir  details  an  improvement  on  that  model  that  accounts 
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rotation  (degrees)  =  -60 
axial  stretch  factor  =  2.  0 


Figine  20.  A  Simulation  of  the  Convergence  Term  of  the  Vorticity  Tendency  Equation  for  the  1826  Del  City  Example 
of  Figure  19a.  The  mesocyclone  orientation  in  this  case  (-60®)  corresponds  to  its  geographic  orientation,  as 
in  Figure  19,  and  not  to  the  orientation  with  regards  to  a  radar. 
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Figure  21. 


rotation  (degrees)  =  -60 
axial  stretch  factor  =  2. 0 


The  Divergence  Field  Contributing  to  the  Convergence  Term  of  Vorticity  Tendency  Shown  in  Figure  20. 
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total  vertical  vorticity 


rotation  (degrees)  =  -60 
axial  stretch  factor  =  2.  0 


Figure  23.  The  Total  Model  Vorticity  Field  From  die  Combinadon  of  a  Uniform  Vorticity  Reid  fw  Rotation  With  the 
Vorticity  Induced  by  the  Non-uniform  Inflow. 
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for  the  observed  elongation  of  many  mesocycione  flows.  This  new  model  does  not  negate  the  concqX  of  a  circular 
mesocyclone.  but  is  rather  an  extension  of  the  original  Ranldne-Combined  vortex  flow. 

Some  exanqries  of  elongated  mesocyclone  flows  observed  with  single-Doppler  radar  are  presented.  These  are 
closely  simulated  by  a  two-dimensional  elliptical  model  that  has  a  simple  kinematics  structure.  The  gmieral  model  is  one 
in  which  vrnticity  and  divergence  are  both  uniform  within  the  core,  a  quality  associated  with  the  amventirmal  “circular'’ 
mesocyclone  model.  One  deviation  from  the  circular  mesocyclone  model  is  a  couplet  of  anticyclonic  voracity  in  the 
potential  flow  regitm  of  a  cyclonically  rotating  elliptical  mesocyclone.  The  intensity  of  the  andcyclonic  couplet 
increases  with  increasing  elongadon  of  the  ellipse. 

An  imptsrtant  caveat  of  ellipbcal  flows  is  that  the  detected  rotadonal  shear  varies  as  a  funcdon  of  the  ellipse 
orientadon  to  the  radar.  The  rado  of  the  maximum  to  minimum  rotadonal  shear  within  the  ellipse  varies  according  to 
twice  the  ellipse  stretch  factor.  For  the  Del  City  exanq>le  discussed  herein,  where  the  major  axis  is  1.8  times  the  minor 
axis  length,  the  detected  radial  shear  would  vary  by  360  percent  around  the  mesocyclone  paimeta-  if  the  flow  were 
perfectly  modeled  by  an  ellipse.  The  problem  is  that  the  established  criteria  for  mesocyclmie  detecdon  assume  that  dte 
mesocycltme  flow  is  approximately  circular.^^  Therefore,  early  detecdon  of  mesocyclonic  drear,  or  detecdon  of  weak 
mesocyclones  at  all,  is  very  much  dependent  on  viewing  angle.  Either  the  critoia  for  nresocyclone  detecdon  need  to  be 
modified,  or  methods  to  infer  peak  shear  in  elliptical  flows  need  to  be  developed. 

Excess  rotadonal  Idnedc  energy  suffers  from  the  same  viewing  angle  jnoblem  as  shear.  Unless  the  eccentricity  is 
known,  the  ERKE  may  be  underestimated  or  ovaesdmated  ctmsideraUy:  by  a  factor  of  2  or  more  for  even  moderate 
eccentricity  widi  a  stretch  factor  of  1.5.  For  a  mid-level,  elliptical  mesocyclone,  where  the  flow  field  consists  solely  of 
rotation,  the  orientation  angle  where  circular  flow  can  be  aviumed  widi  minimal  enor  in  estimating  ERKE  is  about  45 
degrees.  For  the  Del  City,  Oklahoma  storm  examined,  with  an  inflow  from  the  SE  and  the  mesocyclone  major  axis 
oriented  NNW  to  SSE,  the  best  viewing  angle  to  determine  ERKE  is  with  the  storm  generally  east  or  west  of  the  radar. 
Qockwise  (counterclockwise)  of  these  directions  ERKE  would  be  overestimated  (underestimated)  if  the  flow  were 
api^oximated  by  a  circle.  It  would  be  wortiiwhile  to  d^ennine  the  general  orientation  of  elliptical  flows. 

Rotational  kinetic  energy  can  be  determined  accurately  regardless  of  ellipse  viewing  angle  for  an  elliptical  flow  by 
assuming  the  detected  peak  velocities  belong  to  a  circular  flow.  Unfortunately,  RKE  has  not  been  found  useful  for 
mesocyclone  classification  or  for  detormining  sevoe  stom  potential. 

The  model  describes  a  two-dimensicmal  flow  in  what  is,  in  reality,  a  three-dimensional  flew  field.  It  is  easier  to 
apply  constraints  to  the  genoally  horizontal  rotatimi  field  than  to  the  divergence  field,  because  wntii  horizontal 
divergence,  mass  continuity  requires  a  direct  transl^on  to  votical  motion.  Constant  horizontal  mass  flux  for  rotation 
appears  to  be  a  valid  and  useful  assuttqKion  for  the  interpretation  of  the  mesocyclone.  In  examples  presented  herein,  the 
application  of  various  divergence  fields  to  a  uniform  vorticity  field  closely  reproduces  the  observed  single-Dq)pler  flow 
field.  A  good  test  of  the  concept  of  constant  horizontal  mass  flux  for  rotation  is  provided  by  a  multiple-Doppler  aiudysis 
of  the  Del  City,  Oklahoma  mesocyclone  by  Btandes.^  His  analysis  generated  the  actual  vorticity  and  divergence  fields, 
which  can  be  used  as  inputs  to  the  model.  With  a  constant  vorticity  in  the  cote  due  to  rotation,  the  non-uniform 
divergence  field  from  the  multiple-Do{^ler  analysis  is  ^lied.  Interestingly,  the  inflow  field  accounts  for  the  obsCTved 
irregularities  of  the  vorticity  field.  The  resulting  single-Doppler  simulation  closely  resembles  the  obs«ved  teal-life 
example.  This  increases  confidence  that  the  model  is  realistic. 

The  limited  perspective  offoed  by  single-Don>ler  radar  can  make  it  difficult  to  interpret  die  mesocyclone  flow.  It 
is  possible  to  retrieve  more  accurate  flow  estimates  by  tqiplying  more  realistic  models.  With  the  ultimate  goal  of 
replying  automatically  such  a  model  as  presented  here,  much  work  needs  to  be  done.  Present  mesocyclone  detection 
algorithms  for  NEXRAD  do  not  account  for  elliptical  flow  in  their  estimations.  The  results  presented  herein  show  that 
this  shortemning  can  have  dire  consequences  in  terms  of  mesocyclone  detection  and  severe  weatfao'  evaluation. 
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